Maximizing the Neel temperature of fermions in a simple-cubic optical lattice 
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For a simple-cubic optical lattice with lattice spacing d, occupied by two species of fermionic atoms of mass 
m that interact repulsively, we ask what conditions maximize the Neel temperature Tjv in the Mott insulating 
phase at density one atom per site, with equal numbers of the two species. This maximum, fesTi max ' = 
0.15 h 2 /(md 2 ), occurs near the edge of the regime where the system is well-approximated by the usual Hubbard 
model. The correction to the Hubbard-model approximation that produces a "direct" ferromagnetic interaction 
between atoms in nearest-neighbor Wannier orbitals is the leading term that limits how high Tjv can be made. 

PACS numbers: 



One of the next notable milestones in the production of new 
strongly-correlated many-body states with ultracold atoms is 
expected to be the antiferromagnetic Neel phase of two hy- 
perfine species of fermionic atoms in an optical lattice |2|] . 
Important progress towards this goal includes the recent real- 
ization of the Mott insulating phase with fermions [3], and the 
demonstration of controllable superexchange interactions in 
an optical lattice, albeit with bosonic atoms J3]. When the 
optical lattice is sufficiently deep and the repulsive s-wave 
interaction between the atoms is sufficiently weak, the Neel 
temperature Tv for the case of one atom per lattice site can 
be estimated by modeling the system as a one-band Hubbard 
model, and one can analyze the possibility of reaching this 
phase by adiabatically ramping up the interactions and the op- 
tical lattice HI S 0, 0, [^] . The most accessible conditions 
for first producing this ordered phase in an experiment will 
most likely be some compromise between the highest Tjv and 
the highest entropy at the transition S'(Tjv). If the parameters 
of the system, namely the intensity Vq of the optical lattice 
and the s-wave scattering length a s , can be tuned in a per- 
fectly adiabatic manner, then to access the Neel phase only 
requires achieving sufficiently low entropy JHS0]- But m 
the more likely event that there is always some "background" 
heating, so things are not perfectly adiabatic, the phase will 
be more accessible when it occurs at higher absolute temper- 
ature. Thus in this paper we study how the Neel tempera- 
ture Tjv depends on the two tunable parameters Vq and a s as 
one leaves the region where the standard Hubbard model is a 
good approximation to this system. Away from the Hubbard 
regime, theoretical studies have suggested that one may be 
able to access a wealth of phases governed by quantum spin 
hamiltonians HHGH- 

According to quantum Monte Carlo simulations [11] of the 
simple-cubic fermionic Hubbard model, for a given nearest- 
neighbor hopping matrix element t the highest fcsTjv = t/3 
occurs at interaction U = 8 1, while for a given U the maximal 
UbTn — U/20 occurs at t = 0.15 C/. Thus to increase Tjv one 
wants to move to larger t, which means a weaker optical lat- 
tice (smaller Vq), and to larger U, which means larger a s . This 
necessarily moves the system away from the regime where 
it is well-approximated by the usual Hubbard model. The 
mapping from the real system to the Hubbard model uses the 



single-atom Wannier orbitals as the basis states II ll bi ll 211 . The 
standard one-band Hubbard model includes only the lowest- 
energy Wannier orbital at each lattice site and only the on-site 
interaction between two atoms of different hyperfine states oc- 
cupying the same Wannier orbital. 

We find that it is the corrections due to including the inter- 
actions between Wannier orbitals on nearest-neighbor lattice 
sites that are the leading effects that stop and reverse the in- 
crease of Tjv as t and U are increased by decreasing Vq and 
increasing a s . In particular, these interactions produce a "di- 
rect" ferromagnetic exchange interaction favoring neighbor- 
ing sites to be occupied by the same species. These ferro- 
magnetic interactions are of the opposite sign from the anti- 
ferromagnetic superexchange interactions that cause the Neel 
ordering, and thus they suppress Tjv. Within the approxima- 
tions that we make (discussed in detail below) the maximal 
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0.03T r occurs near Vq = 3E r and a s = 0.15 d, 
(irh) 2 1 (2md 2 ) is the recoil energy and d is 
the lattice spacing. For example, E r = 1.4 for 6 Li with 
d = 532 nm, which puts the maximum Neel temperature near 
40 nK, which seems to be well within the reach of current 
experimental cooling techniques. 

This regime of large repulsive a s is attained by approaching 
a Feshbach resonance from the repulsive side. But the atoms 
must scatter repulsively without "falling" in to the weakly- 
bound molecular state. Ref. [31] studied the Mott insulator 
with 40 K at a s = 0.08 d and did not mention any problem 
with excessive molecule formation. It is not clear whether 
this can be increased to the a s = 0.15 d that maximizes Tjv 



1 1_3|] . It is also not clear whether the optical lattice increases or 
decreases molecule formation. The lattice breaks momentum 
conservation, thus possibly opening up channels for molecule 
formation, while in the Mott insulator the atoms are kept apart 
in different wells of the optical lattice, which, naively, reduces 
the opportunities for molecule formation. 

The system we consider is made up of fermionic atoms in 
a simple-cubic optical lattice with a single-atom potential of 
the standard form yfl : 
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This is a separable potential, so the energy eigenstates of a 
single atom in this potential can be chosen to be the product 
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of one-dimensional (ID) eigenstates along each direction. We 
solve for these ID bands and thus obtain the properly nor- 
malized wavefunctions w n (x) of the maximally-localized ID 
Wannier orbitals for each band n 01211 . The 3D Wannier or- 
bitals are then 4>n x n y n z V, z) = w nx (x)w rly (y)w nz (z). 

We will be focusing on the Mott insulating regime with 
density exactly one atom per lattice site and low temperature, 
where the atoms are primarily in the lowest band, n x = n y = 
n z = 0. The nearest-neighbor hopping matrix element to in 
this band is a strongly decreasing function of Vo, the lattice 

strength, behaving as t « 4 7 T- 1 / 2 Er /4 V 3/4 e- 2 ^ /v/Er for 
large Vb [5]. The wavefunction wq(x) of the lowest Wannier 
orbital at a given lattice site is positive and has its maximum 
amplitude at x = 0, the center of the well of the optical lattice 
at that site, while its amplitude is negative and of smaller mag- 
nitude in the nearest-neighbor wells of the lattice. The ratio 
of these amplitudes is one small parameter that is important in 
the approximations we use below. 

We are interested in the antiferromagnetic Mott insulating 
phase in the regime where the on-site repulsive interaction U 
is stronger than the hopping t. We do not treat the limit of 
weak interaction, where the system is a paramagnetic Fermi 
liquid. We assume the atoms equally populate two different 
hyperfine states; as is standard, we will call these two states 
"up" and "down" and treat them as the two states of a spin- 1/2 
degree of freedom. The s-wave repulsive interaction is only 
between atoms of opposite spin. To lowest order in (a s /d), 
we approximate this 2-atom interaction as the standard regu- 
larized contact potential 13 
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FIG. 1 : Approximate phase diagram for filling one fermion per lat- 
tice site. The line marked |J/| = J a is our approximation to the 
ground-state phase boundary separating the antiferromagnetic phase 
at smaller a s from the ferromagnetic phase at larger a s . The ferro- 
magnetic phase is mostly a fully-polarized band insulator, but there is 
a small sliver of polarized Fermi liquid at small Vo between the lines 
marked \Jf\ — J s and "FM to FT'. The line marked "Optimal Tjv" 
indicates where Tn as a function of Vo is maximized for each given 
a s . The dot on that line is our estimate of the parameters that produce 
the overall maximum of Tn/E t , and at that point J/ ~ — J a /4 (see 
text). The U = 14t line is near where the entropy is maximized at 
Tn |5|1 and Tn on this line is maximized at the dot. The to = ti 
line signals when the interaction correction to the hopping becomes 
strong. There is presumably also a paramagnetic Fermi liquid ground 
state in the lower left corner of this phase diagram, but our approxi- 
mations are not well-suited to estimating where this phase is. 
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where r distance between the two particles. The expectation 
value of this interaction energy for two atoms of opposite spin 
occupying the same lowest Wannier orbital is our first esti- 
mate of the strength of the on-site interaction Uon^riii in the 
corresponding one-band Hubbard model: 
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In the Hubbard model, when adjacent sites i and j are each 
singly-occupied by atoms with the same spin, then the hop- 
ping between those two sites is Pauli-blocked. When these ad- 
jacent sites are each singly-occupied by opposite spins, then 
virtual hopping between these sites, treated in second-order 
perturbation theory, allows them to lower their energy and 
thus generates an antiferromagnetic superexchange interac- 
tion J s (Si ■ Sj - \) with J s = At 2 /U. 

The leading corrections to the Hubbard model approxima- 
tion to this system in the regime we are interested in are due 
to the interactions between atoms of opposite spin occupy- 
ing lowest Wannier orbitals on nearest-neighbor sites i and 
j. There are 2 contributions: First, and apparently most im- 
portant in limiting how large Tjy can be made, is the "direct" 



interaction [4] 
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between atoms of opposite spin in adjacent orbitals. This term 
is due to the overlap of the probability distributions of adja- 
cent Wannier orbitals. It raises the energy of the Neel state. 
It thus produces a direct ferromagnetic exchange interaction 
Jf(Si ■ Sj — j) with Jf — —2U nn < that partially can- 
cels the antiferromagnetic superexchange J s that occurs in the 
Hubbard model. It is primarily this ferromagnetic interaction 
that stops and reverses the increase in Tjy as one moves to- 
wards stronger interaction and a weaker lattice while staying 
near the optimal values of U/t. At the global maximum of 
Tn, indicated in Fig. [U we find Jf ~ — J s /4. 

For large enough a s this direct ferromagnetic exchange is 
stronger than the superexchange and thus we have a ground- 
state phase transition to a ferromagnetic phase, as indicated in 
Fig. 1, and discussed more below. [Very near this |J/| = J s 
line, effects due to weaker further-neighbor interaction might 
produce some other magnetically-ordered phases.] 

The "direct" nearest-neighbor interaction (4) also reduces 
the effective U that enters in the superexchange interaction, 
so at this level of approximation our simple-cubic Hubbard 
model has interaction U = Uq — 6U nn , since it is the change 
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FIG. 2: Our estimates of the optimal Neel temperature, Tat, as a 
function of a s /d. For each value of a s , Tjv is maximized by varying 
the lattice depth Vo- We also plot Tjv at the line U — 14t, which is 
near where the critical entropy is maximized (H- 



of the interaction energy due to moving the atom that enters 
in the energy denominator in the superexchange process. 

Also, the interaction generates an additional hopping term 
of the same sign as to: 1511 



ti 



m 



dxw^x)] 2 / dyw%(y)w Q (y + d) (5) 



that is operative when the two sites are each singly-occupied 
by opposite spins. The resulting effective hopping that enters 
in the superexchange process at this level of approximation is 
thus t — to + tj. 

Thus once we include these leading effects due to the 
nearest-neighbor interaction, the effective Hamiltonian in the 
vicinity of the ground state of the half-filled Mott insulator 
has hopping t = to + tj, an effective on-site interaction 
U = Uq — 6U nn , and an additional ferromagnetic nearest- 
neighbor exchange interaction Jf = —2U nn when both sites 
are singly-occupied. To estimate the Neel temperature of our 
system we propose the following approximation: For the Hub- 
bard model without Jf, we have estimates of its Neel tempera- 

ture Tjy (t, U) from quantum Monte Carlo simulations 111 111 . 
This Neel ordering is due to the antiferromagnetic superex- 
change interaction J s = 4t 2 /U between neighboring singly- 
occupied sites. When we include J/ < this reduces this 
magnetic interaction, and we will approximate the resulting 
reduction of Tjv as being simply in proportion to the reduc- 
tion of the total nearest-neighbor exchange interaction: 
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In Fig. 1 we show the lattice strength Vo/E r that maxi- 
mizes this approximation to Tjv for each value of the interac- 
tion a s /d. The highest Tjv occurs near a s /d = 0.15, but the 
system at this value of a s is may be too close to the Feshbach 



resonance and thus not stable against formation of molecules. 
The highest achievable Tjv thus may be somewhere along this 
line at a lower value of a s and thus a stronger lattice Vq. We 
note that a recent experiment y|] has studied a s /d ~ 0.08 for 
40 K, albeit at a temperature well above Tjv, without noting 
any strong instability towards molecule formation. We also 
show on Fig. 1 the line along which U = 14i, since this 
is near where the critical entropy 5 (Tjv) is maximal ]5[], so if 
the system can be adjusted adiabatically this is where the Neel 
phase is most accessible. 

In Fig. 2 we show fcsTjv /E r as a function of a s /d at the 
value of Vq that maximizes our estimate of Tjv, as well as 
at the value of Vo that gives U — 14t and thus is near the 
maximum of 5 (Tjv). Note that in Fig. 2 the horizontal scale 
for a s /d is logarithmic, so Tjv drops rather weakly as a s is 
reduced, meaning that the possible limitation in how large a s 
can be made will not "cost" a lot in terms of the resulting 
reduction of Tjv. 

The approximations we are making are clearly beginning to 
break down in the vicinity of the parameter values that maxi- 
mize Tjv- Thus, although we expect that these approximations 
give reasonably reliable rough estimates of the maximal val- 
ues of Tjv, there are many higher-order effects that we are 
ignoring that may alter these estimates by a little (our calcu- 
lations suggest on the ~ 10% level). At the maximum of Tjv, 
| Jf | is about 25% of J s , The correction to J s due to tj is 
also of roughly this size, but its dependence on a s is much 
weaker, which is why Jf is the important actor in causing the 
maximum in Tjv. 
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FIG. 3: The strongest higher-order process contributing to the energy 
of the antiferromagnetic Mott insulator at the maxima of Tm shown 
in Fig. [U It consists of (1) a nearest-neighbor hop in the lowest 
(S) band, (2) an on-site "pair hopping" of both fermions up to the 
next (P) band, (3) on-site pair hopping back to the S band, and (4) 
a nearest-neighbor hop back to the original configuration. At both 
maxima of Tjv, this process corrects J 3 by about 10%. 

The approximations we have used are those appropriate for 
the Mott insulator, and are based on the inequalities on en- 
ergy scales £rj > U > t, where £o 1S the expectation value 
of the single-particle energy in a lowest Wannier orbital. We 
have analyzed in perturbation theory many corrections be- 
yond those included above. We find that at the maximum 
of Tjv (both the global maximum and the maximum along 
the U — lit line) the strongest next correction is the fourth- 
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TABLE I: The values of the various energies at the two Tn maxima. 
The " J s correction" corresponds to the process detailed in fig. [3] 





Global maximum 

of Tjv 


Maximum 
of T N at U = 14t 


eo (-Br) 


4.2 


5.5 


U (Er) 


0.9 


1.3 


to (Er) 


0.11 


0.07 


tl (E r ) 


0.02 


0.02 


Js (E r ) 


0.08 


0.03 


J J (Er) 


-0.02 


-0.008 


J s correction (-Er) 


0.007 


0.004 



order process illustrated in Fig. [3j it alters J s by about 10%. 
Since our perturbatively-based approximations are breaking 
down near this regime of interest where Tjv is maximized, it 
would be nice to have a more systematic approach that can ob- 
tain more precise and reliable estimates of the phase diagram 
in this regime. For example, quantum Monte Carlo simula- 
tions might be possible for temperatures near TV, although of 
course the famous fermionic "minus signs" may prevent this 
from being feasible in the near term. 

The ferromagnetic phase of this model at strong repulsion 
is mostly a band insulator, with a band gap between the spin- 
polarized bands. But in the weaker lattice regime there should 
also be a partially-polarized Fermi liquid ground state near 
the phase boundary to the Neel state. The transition from the 
fully-polarized band insulator to the partially-polarized fer- 
romagnet occurs when the spin-polarized bands overlap, so 
the system can lower its energy by flipping spins. A sin- 
gle spin flip makes a hole and a doubly-occupied site that 
are each moving freely within the fully -polarized background 
state. At the level of approximation we have used in this 
paper, the hole moves freely with hopping to, so its lowest 
energy is — 6io- The doubly-occupied state costs interaction 
energy Uq + 6U nn and moves freely with effective hopping 
ti = to + 2ti because its motion is the hopping of the flipped 
spin between sites that are both occupied by undipped spins. 
The total energy of this particle-hole pair can be negative 
when Uq + 6U nn < 12(£o + ti) = 6(t + £2); this occurs be- 
low the line indicated in Fig. 1 as "FM to FI" (ferromagnetic 
metal to ferromagnetic insulator). We show this for complete- 
ness, although these ferromagnetic phases at high a s are very 
likely to be inaccessible in experiments with cold fermionic 
atoms in optical lattices. Also, the present approximations are 
probably not very reliable in this regime of large a s /d. 

There is also a paramagnetic Fermi liquid phase at weak 
enough lattice and at weak enough interaction, as well as pos- 
sibly an antiferromagnetic Fermi liquid near it. These phases 
occur well away from the regimes we have focussed on here, 
and the present approximations are not well suited to estimat- 
ing the location of the corresponding phase transitions, so we 
leave that part of the phase diagram as "terra incognita" for 
now. The quantum phase transition between the Neel state and 
the paramagnet should occur in parameter regimes that are ac- 



cessible to the experiments, although it may not be possible to 
see its effects at accessible temperatures, since Tjv decreases 
strongly as this regime is approached. 

Conclusion: We have shown that to maximize the Neel 
temperature one must leave the region of parameter space 
where the Hubbard model approximation for fermionic atoms 
in an optical lattice is well-controlled. We have found that the 
nearest-neighbor direct ferromagnetic exchange is the most 
important correction to the Hubbard model that limits the 
maximal Tn . There are also higher-order corrections to the 
Hubbard model: virtual hopping into higher bands and other 
higher-order processes. For the parameters that maximize Tjy, 
these higher-order terms are smaller than the nearest-neighbor 
terms we include, although not by a large margin of "safety". 

The relative contribution of the higher-order corrections in 
the vicinity of the optimal T/v drops exponentially as one goes 
to smaller interaction a s and thus a stronger optical lattice. 
Thus our results are accurate in the large Vq (strong lattice) 
limit, and should qualitatively capture the phase diagram for 
weaker lattices. For quantitatively more accurate results in 
the weak lattice regime, one needs to resort to more system- 
atic quantum calculations. Of course this is a system of many 
fermions, so it is not clear whether this weak lattice regime 
can be accurately treated in some form of quantum Monte 
Carlo simulations. 
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